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It is shown phenomenologically that the fractional derivative ^ — D°'u of order a of a multifractal 
function has a power- law tail oc |C|~^* in its cumulative probability, for a suitable range of q's. 
The exponent is determined by the condition C,p^ = api, , where ^p is the exponent of the structure 
function of order p. A detailed study is made for the case of random multiplicative processes (Benzi 
ei al. 1993 Physica D 65: 352) which are amenable to both theory and numerical simulations. Large 
deviations theory provides a concrete criterion, which involves the departure from straightness of 
the (p graph, for the presence of power-law tails when there is only a limited range over which the 
data possess scaling properties (e.g. because of the presence of a viscous cutoff). The method is 
also applied to wind tunnel data and financial data. 
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I. INTRODUCTION 

Multifractality for functions was introduced by Parisi 
and Frisch |^ to interpret experimental results of Ansel- 
met et al. on fully developed turbulence g]. We show 
in this paper that multifractality is connected to the tail 
behavior of fractional derivatives (to be defined precisely 
later). First, in Section Owe give a simple phenomeno- 
logical argument predicting power-law tails in probabil- 
ity distributions of fractional derivatives and suggesting 
practical constraints on their observability. In Section III 
we consider the case of a class of synthetic multifractal 
functions, the random multiplicative processes of Benzi 
et al. 11], for which theoretical results are obtained by 
relating them to the theory of random linear maps M. 
In particular, using large deviations theory, we can pre- 
dict the number of generations in a given random mul- 
tiplicative process needed to observe power-laws tails in 
probabilities of fractional derivatives. Numerical experi- 



ments with such processes are presented in Section III C 



In Section IV we analyze high-Reynolds number wind 
tunnel turbulence data and financial data; we show that 
the former, having only weak deviations from self-similar 
behavior a la Kolmogorov 1941, are unlikely to display 
power-law tails at accessible Reynolds numbers. Sec- 
tion ^ gives the conclusions. In this paper we consider 
exclusively multifractal functions; the case of multifractal 
measures, arising, e.g., from chaotic dynamical systems 
or the dissipation of turbulent flows (see, e.g., Ref [^), 
will be addressed elsewhere. 



h < /imax < 1) there is a continuum of sets Sh of fractal 
dimension D{h) such that, when x belongs to Sh, the 
function u has a Holder-like behavior with exponent < 
h < 1 m. In terms of increments this is expressed as 



\Su{x, r)\ ex \r\ 



0, X e Sh, 



(1) 



where Su{x,r) = u{x + r) — u{x). If we take D{h) 
to be the covering dimension (see Ref. ||] Section 8.5.1 
and references therein) and use instead the codimension 
F{h) = d — D{h), where d is the dimension of space, 
and then "thicken" the set Sh by covering it with balls 
of small radius p, into the set S^^, we have 



Prob{x e Sj"^} ex p 



F(h) 



p^Q. 



(2) 



For the case of an ergodic random function, this proba- 
bility can also be interpreted as a fraction of space. Note 
that by using codimensions the extension from the one- 
dimensional to the multi-dimensional case is straightfor- 
ward (just replace scalars used throughout this paper 
by vectors). For more background on multifractals see 
Ref. [^ and references therein. 

We turn now to fractional derivatives. The definition 
used here for the fractional derivative D'^u{x) of order 
< a < 1 of a function u{x) is to multiply its Fourier 
transform u{k) by |fc|". This might be called, more cor- 
rectly, the fractional (negative) Laplacian of order a/2. 
In physical space our fractional derivative may be writ- 
ten, in the one-dimensional case, as 



D"'u{x) 



^^sin^r(a)P.V. / 

TT 2 ^ ^ J 



u{x + r) — u{x) 



■\a+l 



dr. 



(3) 



II. PHENOMENOLOGY 



A (homogeneous random) function u{x) with a; G IR 
is called multifractal if (for some real interval < /imin < 



where T{a) is the Gamma function and P.V. denotes a 
Cauchy principal value. The Fourier-space definition is 
very convenient for periodic functions. When working 
with non-periodic experimental data, windowing may be 
necessary before the Fourier transformation is applied. 



Note that other definitions of the fractional derivative, 
for example that based on the ordinary derivative of the 
the Riemann-Liouville fractional integral , do not pro- 
duce substantially different results, as far as our work is 
concerned. 

It is clear from (m that, at a point xq at which u{x) has 
Holder exponent h > a the fractional derivative D°'u{xq) 
is finite (provided the function u is, say, bounded). If a > 
h the fractional derivative will generally be infinite. This 
is also the case with the definition of fractional derivatives 
used in Ref. B . But what takes place in the neighborhood 
of such a point? Let us first consider the case of an 
isolated non-oscillatory singularity of exponent h at x. 
That is, we have \u{x + r) — u{x)\ oc \r\'^ for small r. 
By substitution into (H) it is easily checked that, with a 
suitable constant B > 



\D''uiy)\ ^ B\y - 



(4) 



a relation which could have been guessed by simple di- 
mensional analysis. Since h—a < 0, the fractional deriva- 
tive becomes very large near x. 

What happens if instead of isolated singularities we 
have multifractal non-oscillatory singularities? By a 
standard argument of multifractal analysis, in the phe- 
nomenological presentation of Ref. [Q (see also Ref. ||] 
Section 8.5), we still use (0) for determining the contri- 
bution of points x G Sh to D^ud/) calculated at points 
which are near Sh without belonging to this set. It is now 
very simple to estimate the probability to have |I?"u| > £, 
for large positive ^. From (y), the contribution from a 
given h < a to \D°'u{y)\ will be greater than ^ provided 
that 



x|<l| 



By (0), this has probability 



Prob {{D^ul > C} oc |y - x^ ^ ' oc ^ 



(5) 



(6) 



This is a power law in ^. As we vary the Holder exponent 
h, the exponent of this power law changes. When we sum 
the contributions of the various ft-'s (with suitable regular 
weights dfj,{h) which we need not specify) we get an inte- 
gral which, when evaluated by Laplace's method, gives us 
leading-order power-law behavior for the tail probability: 



Pioh {\D°'u\ >^}(xC^*, 
F{h) 



p^{a) = inf 



h<a a — h 



(7) 
(8) 



We shall now show that p^, can be easily found from the 
scaling exponents Cp of structure functions. We remind 
the reader that, when u{x) is a homogeneous multifractal 
function, the moments of its increments, called structure 
functions, are defined as 



S,{r)^{{\Su{r)\r), p>0. 
They follow power laws at small r's: 
Sp{r) ex \r\i-, 



(9) 



(10) 



where the up convex function ^p and the down convex 
function F{h) are Legendre-Fenchel transforms of each 
other i § 

F{h) = snp{Cp-ph), Cp = mi{ph + F{h)). (11) 



Using(ni|) in (H) and interchanging the orders of the in- 
fimum and the supremum, we obtain 



. f Cp~ph 

Pi, — sup mt -^ — . 

p>0 h<a a ~ h 



(12) 



The infimum over h < a is — cx) if (^p—pa < 0, it equals p if 
(^p—pa = and is less than this value if (^p—pa > 0. Thus 
Pi,{a) is the supremum of all p's such that (p — pa > 0; 
in other words, it is the solution of 



(p. =P*(a)a, 



(13) 



which is unique in view of the convexity of ^p and vanish- 
ing for p = (a consequence of (0)). Thus, the exponent 
p*{a) of the tail probability of the fractional derivative 
of order a can be obtained from the graph of (p by the 
simple construction shown in Fig. n. a line through the 
origin of slope a intersects the graph at p = Pi,{a). 




FIG. 1. Geometrical determination of the exponent p*(q) 
of the power- law tail from the graph of (^p, the struc- 
ture-function exponents, /imax is the maximum Holder ex- 
ponent; hmin, the minimum Holder exponent, is the (alge- 
braically) smallest slope, not shown. 

This construction has a number of interesting conse- 
quences which are now presented. First, recall that 



dp 



(14) 



p=p* 



(This follows from the observation that F{h), if down 
convex, can be recovered from C.p by another Legendre 



transformation.) If there is a minimum Holder exponent 
/imin and maximum exponent /imax, it is clear that a has 
to be between those two values for p^, to exist. If a < /imin 
we do not expect any power-law tail at all (at least not by 
the mechanism considered here) . If a > /imax we expect 
D°'u{x) to be infinite almost everywhere. D°'u{x) can 
then be made finite if we assume that there is some ul- 
traviolet cutoff, e.g. a viscous cutoff. If so, the tail prob- 
ability of D°'u{x) can be estimated by the same kind of 
phenomenological arguments used for the tail probability 
of velocity gradients in turbulence P,p|pX|] . 

Next, we can use the construction to estimate how far 
the power-law tail (if present at all) is expected to ex- 
tend when multiscaling holds only over a finite range of 
scales. By (Q) the fractional derivative stemming from 
the Holder exponent h^, is given by 



\D^u{y)\ 
C{a 



a 



P* dp 



(15) 
(16) 



C{a) will be called the multifractality parameter. In the 
next sections we shall see that the presence of a power- 
law tail requires in practice 



nC{a) > 10, 



(17) 



where n is the number of octaves in the scaling range 
(e.g. in the inertial range for turbulence). If the range 
of spatial scales is somewhat limited, it is crucial that 
C{a) be as large as possible. In those instances where 
the graph of Cp is close to a straight line through the 
origin, the two terms on the r.h.s. of ( |l6[ ) will nearly 
cancel and C(a) will be very small. Clearly, having a 
graph departing strongly from a straight line will help in 
seeing power-law tails. 

All the material presented in this phenomenological 
section is, by definition, "soft" . Is it possible to give 
harder evidence? One way is to work with Burgers tur- 
bulence in the limit of vanishing viscosity (Ref. ||ll| and 
references therein) . The singularities (mostly shocks) are 
then isolated and it is easy to make our phenomenolog- 
ical arguments rigorous. Solutions to the Burgers equa- 
tion are however bifractal, not truly multifractal. In the 
next section we discuss a simple example of a truly mul- 
tifractal function. 



differs somewhat from that of Ref. g] but is basically 
equivalent. 

We first define the local rmp as a real random function 
on the real line. Let 4>{x) be an indefinitely differentiable 
function with rapid decrease at infinity and vanishing 
space integral. In practice we shall take the Mexican hat 
function of width a used in Ref. IHl 



da; 2 



exp 



X 

2^ 



(18) 



Let ?Ti be a random variable, called the random multi- 
plier, with symmetric distribution and finite moments of 
all orders, characterized by its cumulative distribution 



P{ri) = Prob {m > 77} 
and its probability density 

pim = — 1 — =p[-m- 

drj 



(19) 



(20) 



We shall also need independent identically distributed 
copies of rn, denoted m„ and m„ [n — 0,1,...). 
(Henceforth the words "identically distributed" are un- 
derstood.) We define now a sequence of random functions 
Wnix) recursively by 

Wo{x)^4>{x), (21) 

w^„+i(2:) = 0(x) + m[l^w';l\2x) -f m'^^^ w':;^\2x - 1), (22) 

71 = 0,1,..., 

where w„ [x) and Wn {x) are independent copies of 
Wnix). The almost sure pointwise limit for ri -^- 00 of 
Wnix), if it exists, is denoted by 'w{x) and satisfies obvi- 
ously 

w{x) ^^<i>{x) + m(i)u;(i)(2x) + m^^'^w^^\2x - 1), (23) 

where = designates equality in law. It may be checked 
that our random map definition is equivalent to that 
given in Ref. H], which involves an infinite series. 

Assuming the limit to exist, we now define the global 
rmp as 



i — + CXD 



l{x) = ^ w'^''^(a; -f Xj), 



(24) 



III. THE CASE OF THE MULTIFRACTAL 
RANDOM MULTIPLICATIVE PROCESS 

In 1993 Benzi et al. introduced an explicit method 
for constructing multifractal random functions with arbi- 
trary (p. Such functions are constructed by iterating suit- 
able random maps and will be called here "random multi- 
plicative processes" (rmp). Our method of construction 



where the w^'^'{x) are independent copies oi'w{x) and the 
Xi are Poisson-distributed points on the real line with 
density c. 

The construction of the global rmp from the local one, 
which differs slightly from that of Ref. |H], guarantees 
statistical homogeneity. Denoting statistical means (ex- 
pectation values) by angular brackets, we can relate the 
characteristic functionals of u{x) and w{x). Indeed, we 
have 



exp 



exp < c 



t(p{x)u{x) dx 



R 



exp 



iip[x)w{x + y) dx\ ) — 1 



R 



dy . 
(25) 



where (p{x) is a test function and jL denotes integration 
from — oo to +00. (To establish (|2^ ) it is convenient to 
assume that there are N points Xi distributed uniformly 
in the interval [— i,+L] and then let N -^ oo, i — > 00 
and N/{2L) ^ c). A special case of ( P5|) relates just the 
characteristic functions (which for u does not depend on 
x): 



be random matrices Q. In its scalar form (^l]) has at- 
tracted the attention of statistical physicists interested in 
disordered systems or finance |lj-[lj] . Results on the ex- 
istence of invariant measures for Kesten maps and more 
general nonlinear random maps are found for example 
in Refs. pa-n^. Here, we are only interested in the 
scalar linear map defined by (31), which we shall call 



the "Kesten map" associated to the rmp. 

We now summarize a few important results for Kesten 
maps without giving proofs, but we shall give occasional 
hints. Let 



A = (In \m\) 



(33) 



(e' 



exp 



((■ 



,izw(y) 



1 dy 



(26) 



The fractional derivative D°'u{x) is obtained by replac- 
ing, in (|2|), w{x) by D'^wix). Applying D" to (|3|), we 
find that D°'w{x) satisfies 



Wa[x) 
Wa{x) 



'= 0„(x) -F 2"m(i)u;W(2a;) -f 2"m(2)u;(^2)(2x - 
= D°'w{x), 4>„{x) = D"(j){x), 



I) 



which is basically the same as (p3) with changed into 
D^cj) and the random multiplier rescaled by a factor 2". 
Various theoretical results for structure functions were 
obtained for rmp's in Ref. 0]. In particular it was shown 
that the structure function of order p > has (with our 
notation) the scaling exponent 



be the Lyapunov exponent of the map. If A > the 
w„'s run away to infinity almost surely. (In \fn1m2 ■ ■ ■ mn\ 
grows approximately as An.) If |m| < a < 1 al- 
most surely, then the series (p^ ) converges to a value 
< 1/(1 — a). Hence, the tail behavior of w at large values 
is trivial. When A < and there is a non-trivial tail, 
the nature of the invariant measure (the distribution of 
(27)u)) depends in particular on the average slope {\m\) of 
(28)the map. If (|m|) > 1 it may be shown that the invari- 
ant measure is absolutely continuous, i.e. that there is 
a finite probability density for w. (This may still hold, 
but only in special instances, for {\m\) < 1.) The most 
important result is that, for (|m|) > 1 and A < 0, there 
is a power-law tail: 



Prob{|u;| > £,} (x C 



e 



Cp = -ln2{\m\P) 



(29) where Pi, is the single positive number such that 



Here, we need new results for probability distributions. It 
turns out that there is a class of space-independent linear 
random maps (Kesten maps) whose probability distribu- 
tions are very closely connected to those of the rmp's. 
They are discussed in the next subsection. 



A. Kesten's random maps 

Consider the sequence of random variables defined re- 
cursively by 



(|mr)=l. 



(34) 



(35) 



Wo 



1, 



Wn+l = 1 +TO„W„, 

n = 0, 1,.. ., 



(30) 
(31) 



(Take the absolute value of (|3l|), raise it to the pth power, 
average and neglect the contribution of the additive 1 
when p is close to p* and just below, so that this moment 
comes predominantly from the tail behavior of w.) 

Now we give an alternative derivation of this power-law 
behavior, based on large deviations, which is in spirit very 
close to the phenomenological derivation of power laws in 
Sectio n [ll[ F or background on large deviations see, e.g., 
Refs. |18| , |19[ (an elementary introduction may be found 
in Section 8.6.4 of Ref. §). 

Consider the nth term on the r.h.s. of (p2). Since the 
m variables have symmetric distributions, it may written 



where the 77i„'s are independent copies of the random 
multiplier m, with the same definition as used above. 
The limit of w„ for n — > 00, if it exists, has the same 
distribution as 



± |toi||to2| 



__ I 9 — (ai+a2 + ...+a„ 



= -ln2 \mi 
(36) 



w 



1 + mi + rni?77,2 + . . . + mim2 ■ ■ ■ to„ -t- 



(32) 



where the plus and minus signs arc taken with equal prob- 
abilities 1/2. 

First, let us assume that just one term dominates the 
tail behavior of w: 



Kesten has obtained a number of important results for 
a more general class of maps where the multipliers may 



w 



±|mi||TO2| 



(37) 



(or, more precisely, that the right scaHng is obtained; we 
shall come back to this later). Large deviations theory 
(expressed here somewhat loosely) tells us that, for large 
n and given /i < 0, 

Prob {ai + aa + • • • + a« « nh} - 2'"^^''^ , (38) 

where F{h) is the Cramer (or rate) function, defined here 
to be positive to ensure consistency with the definition 
of the codimension used in Section O (this is also why 
we use powers of 2 and the minus sign in the definition 

of fli). 

From (g^ and ( p8|) , and after integration over all pos- 
sible /I's, we obtain 



Prob{|w| >C}(xr^% 
F{h) 



inf 

/i<0 



Similarly, for p > 0, we find 

((|mi||m2|---|m„|)P)-2-'<-, 



Cp = - ln2 



= inf [ph 

h 



F{h)). 



(39) 
(40) 

(41) 
(42) 



Note that ( [i^ ) is just Cramer's result that F{h) is the 
Legendre transform of the logarithm of the characteristic 
function of the variables a^. Note also that (p9|), (40), 
(IH) and (H) resemble (0), (|), (|) and (^ estabhshed in 
Section OTlt is easily shown by manipulations resembling 
those of Section O that Pi, has exactly the value predicted 
by (pq) and that the minimum in (Bfl) is achieved at 



K 



(|?77,|^* ln2 |r7T,|) 



(43) 



Obviously, in order for p^ to exist, negative values for 
h must be permitted or, in other words, the variable |?7i| 
must be allowed to take values in excess of unity. 

In evaluating the tail of the series (JS^) we have so far 
assumed that it suffices to take a single term of suitable 
order n. Two questions arise: (i) how large should n be 
and (ii) what happens if we take the whole series? The 
answer to the former follows from (p6|), ( p7|) and (38): 
since \w\ « 2""'**, we need to have n of the order of 
ln2 |w|/|/i^|. We observe that, the larger \hi,\, the faster 
(in n) large |w|-values will be reached. Consistently with 
the definition given in Section |l| we call C = |/i*| the 
multifractality parameter. For the latter question we just 
briefly report in words our findings which are a bit too 
technical for the present paper: each of the terms in (p3) 
is, except for a random sign factor, an exponential of a 
sum of Gi terms. These sums perform a random walk. 
If, for large n, we put the condition that there should 
be a large deviation such that the slope of the walk is 
approximately h ^ (a) , then the walk is itself close to a 
straight line of slope h (with small Gaussian fluctuations 
around it). One then finds the same scaling result as 
given by (p3) and it may be checked that the inclusion 
of more than one term just affects the constant in front 
of the power law. 



B. Multifractal properties of the random 
multiplicative process 

The fact that the rmp's, as defined at the beginning 
of Section III, are multifractal has already been proved 
in Ref. [||. The (p function given there (our eq. (p^)) is 
precisely the same as (^2|) for the Kesten map. This is, 
of course, not accidental since the rmp was introduced 
by analogy with random multiplicative cascade models, 
for which a simple relation is known to exist between 
their multifractal properties and the large deviations of 
the logarithms of the multipliers (see, e.g.. Section 8.6 of 
Ref. §). 

Now we turn to tails of the fractional derivatives of 
rmp's and show that they have precisely the properties 
conjectured in Section O. We shall limit ourselves to 
instances where the mean slope of the associated Kesten 
map is greater than one or, equivalently using (pa), p+ < 
1. Let us make a preliminary remark regarding rmp's for 
which the largest multiplier is exactly unity, as is the case 
of the example considered in Section 4 of Ref. g] , when 
account is made for a factor v2 in equation (9) of that 
reference. In such a case, neither the associated Kesten 
map nor the rmp itself can have a power-law tail for w: 
if the multiplier 1 has finite probability, say p < 1, after 
n iteration, w is at most 0(n) with probability p". If 
we now take a fractional derivative of order a > 0, we 
have already noted that this rescales the multipliers by a 
factor 2" so that power-law tails are no more ruled out. 
In the Kesten map associated to D^w there will be a 
power-law tail with exponent —p^, given, in view of (|3£ 
by 



{2P* 



= 1. 



(44) 



By (B3) this is equivalent to Cp* = P*ct, which is identical 
to (p^. For the full rmp we do not know a rigorous 
derivation of this relation but we shall now show that if 
D°'w{x) has, for all x, a power-law tail probability with 
exponent — p* (0 < p* < 1), then Pi, is given by ( [44| ) and 
the same exponent applies to the global homogeneous 
process D'^u. 

For this, following rather closely an argument given in 
Ref. fECi, consider 



Ka{z,x) 



^■izD^w{x) 



(45) 



the characteristic function of the fractional derivative of 
the local rmp. By ( |27| ) and the independence of wS^\ 
mP'\ w^^^ and w'^^-', the following integral equation is 
obtained 

K^{z,x) = (e^D-^i^.)'^ j K^^2"y,2x)piy)dy 



R 



X / K^{2''y,2x-l)p{y)dy, (46) 



R 



where p{-) is the probabihty density of in. By the as- 
sumption made, there is a small-z expansion of the form: 



Kaiz,x) = l + \z\P*gaix)+h.o.t. 



(47) 



where h.o.t. (higher order terms) stands for o(|z|''* ). Ex- 
panding ([iq), collecting all O i\z\P*) terms and integrat- 
ing overs X, we find that 



9a^9a / 2P*"\y\P*p{y)dy, 



(48) 



R 



where ga = J^ ga{x) dx. Eq. ( [44|) follows immediately 
from (Bq). From a small-z expansion of (Eq), with u and 
w replaced by D^u and D'^w, respectively, it follows that 
D'^u has the same tail exponent — p*. 

Finally, we briefly indicate the changes in tail proba- 
bilities which occur when the order of fractional differ- 
entiation a varies. For simplicity we assume that m is 
bounded and, without loss of generality, that the largest 
value rrimax = 1 (otherwise change a into a — ln2 mmax)- 
For a = we already noted that no power-law tail 
can be present. The down convexity of the function 
Pi, 1-^ 2"P* (|to|p*) implies that (H) has at most one so- 
lution other than p* = 0. Whether or not is has one 
depends on the Lyapunov exponent of the Kesten map 
with multiplier 2"to. If < a < — (ln2 |m|) there is a 
power-law tail. If a > — (ln2 \m\) there is runaway. 



C. Designing numerical experiments with the 
random multiplicative process 



The general class studied will be 

_ J ±mi, each with probability g/2; , . 

|_ ±1, each with probability (1 — q)/2, 

where < mi < 1 and Q < q < 1. It is easily shown 
that for the Kesten map associated to u the Cp's and the 
multifractality parameter are given by 

Cp--ln2(gmf + l-g), (51) 

ln2 {qml' + 1 — q) qm^* hi2 mi 



C{a) 



P* 






1-9 



(52) 



where a and Pi, are related by (M). From ( |5l| ) we can 
calculate the minimum and maximum Holder exponents: 



A 



10-' 



0, 



-grln2mi. 



(53) 
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FIG. 2. The (local) Roman rmp [H] and its (^p function (as 
inset). 

In this Section we perform numerical experiments with 
a class of rmp's, called dyadic, where the multipliers take 
only two values (and their opposites). The example con- 
sidered in Section 4 of Rcf. H, here called the Roman 
rmp, is an instance. It has 



_ / ±2-^/^, each with probability 7/16; 
±1, each with probability 1/16. 




(49) 



FIG. 3. (a): distribution of fractional derivative of order 
a — 0.2895 (p* = 0.5) for the Roman rmp [y]. Number of 
iterations (levels) n in the Monte-Carlo simulation with lO'^ 
realizations as labeled. The slope of the line is the predicted 
asymptotic behavior. No power-law is observed, (b): distri- 
bution for the associated Kesten map ( |3l| ) with 10 realiza- 
tions. Power-law behavior is observed only when the number 
of iterations n > 5 x 10^ . 



The simulations reported hereafter have n < 16 itera- 
tions of the rnip (but much more for the Kesten maps). 
The width of the Mexican hat is always cr = 0.1. For the 
construction of the local rmp, the simulation interval, 
from —0.7 to 1.3, is chosen such that the local rmp w{x) 
is down to (single-precision) roundoff level at the edges 
of this interval. The resolution is 2"*^"+'*^ (the -1-4 is 
needed because the width of the Mexican hat is substan- 
tially smaller than unity) . The global rmp is constructed 
by adding typically one hundred independent Poisson- 
shifted local realizations with a density of the Poisson 
distribution c — 0.5. Fractional derivatives are calculated 
as multiplications by \k\" of the (discrete) Fourier trans- 
form (windowing not needed because of the rapid falloff 
at the edges). Cumulative probabilities are calculated 
from the global rmp using rank ordering to avoid bin- 
ning pip^ ]. Cumulative probabilities for Kesten maps 
are obtained by Monte-Carlo simulation of the random 
map with a number of iterations ranging from n = 13 to 
n = 10^ and typically lO** realizations. 



w{x) 




FIG. 4. (a): the (local) optimal rmp defined by (M) and 
its ("p function (as inset).; (b): fractional derivative of order 
a = 1.368 applied to this process; note the spikes. 

Fig. g shows a graph of one realization of the (local) 



Roman rmp for n — 13. This has /imax = 7/24 Ri 0.29. 
Fig. ||(a) shows the cumulative probability of the frac- 
tional derivative for an a chosen to give p^ = 1/2 and 
n = 16. No power-law tail is observed. The associated 
Kesten map does reveal the power-law tail but only after 
about 5000 iterations (in the rmp this would correspond 
to the smallest scale being 2"^°"°!). The reason for this 
slow convergence is the exceedingly small value of the 
multifractality parameter C — 0.0022991. Hence, the 
Roman rmp is not a good candidate for seeing power-law 
tails in fractional derivatives. 




FIG. 5. As in Fig. 0, but for the optimal rmp (M) and 
a = 1.368 (p* = 0.876). Note the conspicuous power-law tail 
after about 13 to 16 iterations. 

For obtaining better candidates we systematically 
search through parameter space, varying mi, q and a 
and trying to maximize the multifractality parameter 
C(a) given by (|5^). In practice, after delineating the ac- 
ceptable ranges of these parameters, we randomly choose 
3 X 10^ triplets (mi,(7, a) and each time calculate C(a). 
We thus obtain about half a dozen triplets for which C[a) 
is close to unity, i.e. more than 400 times larger than the 
value for the Roman rmp. We then select the one for 
which the Kesten map displays a conspicuous power-law 
tail for n as low as possible. This turns out to be n « 13 
(around n « 10 a power-law tail begins to emerge). We 
call this the "optimal rmp" . Its parameters are: 



TOi = 0.04, q = 0.6, a = 1.368, C{a) = 0.9867. 

(54) 

Fig. Q shows the local optimal rmp and its fractional 
derivative of order a = 1.368, which gives p* = 0.876. 
Note that the graph of (p (inset on Fig. §(a)) is much 
more curved than that for the Roman rmp. Fig. |^(a) 
shows the cumulative probability which is seen to dis- 
play a power-law tail already for n = 13 and n = 16. 
Note that the cumulative probability for the associated 
Kesten map Fig. ra(b) is a somewhat wiggly power law 
due to the presence of log-periodic subdominant correc- 
tions |l^] ; this phenomenon is less pronounced if we use 
lower values of p* (e.g. 0.5, not shown). 

IV. ANALYSIS OF EXPERIMENTAL DATA 
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FIG. 6. Distribution of fractional derivative of order 
a = 1/3 for the Modane turbulence data. The power law 
with exponent — p* = —3 is the predicted tail behavior at 
infinite Reynolds number. In fact, no scaling region is ob- 
served. Inset: corresponding (^p graph which is close to its 
Kolmogorov 1941 value, a straight line of slope 1/3. 

As already stated in the Introduction, multifractal- 
ity was introduced to interpret experimental turbulence 
data, specifically the up convex bending of the graph of 
(p. Some of the best high- Reynolds turbulence data have 
been collected at the ONERA wind tunnel by the Greno- 
ble group from LEGI. The Taylor-scale based Reynolds 
number is Rx ~ 2700 with about three decades of scaling 
(for definitions see, e.g., Ref. [p|). We use (longitudinal) 
velocity data obtained from this group to check for the 
possible presence of power-law tails in fractional deriva- 
tives. A total of 2.3 X 10^ data points (in 12 segments of 
variable lengths) is analyzed. To calculate the fractional 
derivatives of these non-periodic data we use Hann win- 
dowing (see, e.g., Section 13.4 of Ref. (23|). The different 
segments are processed separately and the union of all 
fractional derivative data is then rank ordered to produce 
cumulative probabilities (using just a single segment does 



not produce significantly different results; as we shall see 
the problem we shall encounter has little to do with noise 
which can be reduced by using more data). 

Fig. ^ shows the probability of the fractional derivative 
of order 1/3, a value chosen because, by Kolmogorov's 
four-fifths law, the corresponding p^, should be exactly 
three 0. Fig. M shows the case a = 13/36, which cor- 
responds to Pi, = 1/2 (when using the lognormal law 
Cp = (p/3) + (m/18)p(3 — p) with ^ = 0.2, an excellent 
approximation for p up to at least six). For both in- 
stances we show the nearly straight graph of (p over the 
relevant range oi p's. In neither case do we see any trace 
of the predicted power law tail. 
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FIG. 7. As in Fig. |6| but for q = 13/36 (p* « 1/2). 

Here, a digression is in order. The fractional derivative 
analysis of this turbulence data was actually made by us 
before the work on random multiplicative processes re- 
ported in Section [II. The latter was motivated in part 
by our desire to interpret the failure to see power laws in 
experimental data. Of course, now it is clear why we do 
not see such power laws. The multifractality parameter 
of the experimental data, using the lognormal approxi- 
mation, is 



C(a) 



M 



6 



(55) 



For a ^ 1/3 we have C « 0.033 while for a = 13/36 
we have C w 0.0055. These values are far too small 
for power-law tails to be observable with a three-decade 
long inertial range. Actually, what we observe in Fig. o 
and at high values of the modulus ^ of the fractional 
derivative is due to the viscous cutoff and can be pre- 
dicted by the argument of Ref. , adapted to fractional 
derivatives. 

Higher Reynolds number data with four decades of 
scaling will soon become available thanks to the use of 
large-scale cryogenic Helium facilities [M. This increase 
is unlikely to be sufficient to see the predicted power laws. 
Alternative methods, based on the processing of the local 
dissipation measure by fractional derivatives or integrals 



should also be explored. 

We now consider financial data for which the obser- 
vation of power-law tails for fractional derivatives could 
be easier, although the data are far scarcer than in tur- 
bulence. Stock market indices are generally believed to 
display multifractal scaling. We refer, for example, to 
http://www.ncrg.aston.ac.uk/~vicenter/financial.html 
which contains a number of references to multifractal (or 
spuriously multifractal) behavior of financial time series. 
We use S & P Futures intraday (minute-by-minute) data 
with linear interpolation for the approximately 12 % of 
missing data points due to lack of activity. The total 
number of points including interpolated data is 1.88 x 10^. 
We analyze their logarithms (because returns are defined 
as differences of logarithms). Hann windowing is used 
again (higher-order windowing makes no substantial dif- 
ference.). Structure functions (not shown) of order up 
to p — 4: display almost three decades of scaling. The 
corresponding ^p function is shown as an inset on Fig. |g. 
For a = 0.348 we show the cumulative probability (by 
rank ordering) in Fig. g. The Cp graph is much more 
curved than that for turbulence data. The multifractal- 
ity parameter, C(0.348) w 0.30, being almost one order 
of magnitude larger than for the turbulence data, it is 
not surprising to observe a power-law range (oc |CI~^'^) in 
the cumulative probability. Let us finally point out that, 
in one respect, such financial data are closer to Burg- 
ers turbulence (which easily displays power-law tails for 
fractional derivatives) than to Navier-Stokes turbulence. 
Indeed, there are occasional strong discontinuities, such 
as financial crashes. 




FIG. 8. Distribution of fractional derivative of order 
a = 0.348 (p* = 3.5) for S & P 500 Futures intraday financial 
data (minute-by-minute, 1982 April 21 - 1999 February 26). 
The bump around ^ = 2 is mainly coming from the vicinity 
of the crash of October 1987. Inset: corresponding (p graph 
and line of slope a. 



V. CONCLUSION 

We have shown in this paper that multifractality of a 
function has a signature in the tail behavior of proba- 
bilities of fractional derivatives. This signature may fre- 
quently be invisible, since power-law tails emerge only 
if the number of scaling octaves times the multifractality 
parameter exceeds a value found empirically to be around 
10. We recall that the multifractality parameter C is a 
measure of how strongly the functions departs from be- 
ing self-similar. Incompressible three-dimensional turbu- 
lence having C ~ 1/30 is not a good candidate for seeing 
such tails. There are many other turbulence-like prob- 
lems displaying multifractality, such as passive or active 
scalars and magnetohydrodynamic turbulence. Most of 
these should have larger C's. 

It is known that the standard way of detecting multi- 
fractality, by structure functions, can produce spurious 
multifractality. One reason is that self-similarity does 
already produce structure functions with scaling, albeit 
with (p linearly proportional to the order p. If there 
is a mechanism, such as contamination by subdominant 
terms [^ or very slow convergence |^,0 , which slightly 
changes the apparent value of Cp, a self-similar function 
may look multifractal. Such a phenomenon is ruled out 
with our new method which detects only very robust 
forms of multifractality. 

We mention here some open mathematical issues. Jaf- 
fard has recently shown that multifractality is a (topo- 
logically) generic phenomenon in many function spaces of 
the Sobolev family such as Besov spaces, irrespective of 
any underlying dynamics ^^ IQ ■ Can such a proof be 
extended to show genericity of power-law tails? For ran- 
dom multiplicative processes it should be possible to give 
a rigorous proof of some of the results presented some- 
what heuristically in Section III. We notice also that 



the functional random map formulation ( |23| ) can be used 
to investigate further the statistics of individual Fourier 
modes, studied in Ref. [Q. 

Finally, returning to three-dimensional turbulence we 
cannot resist asking: since experimentally this turbulence 
does not "test positive" when looking for power-laws tails 
in fractional derivatives does it have other features which 
would reveal themselves only at vastly larger - and there- 
fore perhaps inaccessible - Reynolds numbers? Could it 
be that the velocity itself (without taking any derivative) 
or that velocity increments have power-law tails in their 
probabilities, as suggested for example in Ref. Q? This 
cannot be ruled out completely. For example, if we ac- 
cept the lognormal model (with /x « 0.2) also for large 
values of p (for which the exponent h is negative), we 
find that (p vanishes for p « 33. This implies a power- 
law tail with exponent —33. The multifractality param- 
eter is C ~ 0.37. Such an algebraic tail, if it exists, 
would require about 30 octaves or 9 decades of scaling. 



that is R\'s of about eighteen milhons. Getting there re- 
quires the crossing of a substantial "turbulence desert" , 
but a small one in comparison with the great high-energy 
desert predicted by Giorgi's and Glashaw's grand unified 
theory ||. 
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